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In Cerium-based heavy électron materials, the 
4/ electron's magnetic moments bind to the 
itinérant quasiparticles to form composite heavy 
quasiparticles at low température (T). The vol- 
ume of the Fermi surface (FS) in the Brillouin 
zone incorporâtes the moments to produce a 
"large FS" due to the Luttinger theorem. When 
the / électrons are localized free moments, a 
"small FS" is induced since it contains only broad 
bands of conduction spd électrons. We have ad- 
dressed theoretically the évolution of the heavy 
fermion FS as a function of T, using a fîrst prin- 
ciples dynamical mean-field theory (DMFT) ap- 
proach combined with density functional theory 
(DFT+DMFT). We focus on the archetypical 
heavy électrons in Ceirins, which is believed to 
be near a quantum critical point. Upon cooling, 
both the quantum oscillation frequencies and cy- 
clotron masses show logarithmic scaling behavior 
ln(To/r)) with différent characteristic tempér- 
atures Tq — 130 and 50 K, respectively. The en- 
largement of the électron FS's at low T is accom- 
panied by topological changes around T = 10 ^ 
50 K. The resistivity cohérence peak observed at 
T ~ 50 K is the resuit of the compétition between 
the binding of incohérent 4/ électrons to the spd 
conduction électrons at Fermi level (Ep) and the 
formation of cohérent 4/ électrons. 

The FS volume has been a sensitive probe of the 
character, localized or itinérant, of the heavy fermion 
System.^ Intensive efforts have been devoted to the study 
of the quantum phase transition leading from a small to 
large FS at strictly zéro température. While the FS, 
as a surface of discontinuity in the momentum distribu- 
tion function, is sharply defined only at zéro température, 
expérimental probes such as the angle-resolved photoe- 
mission spectra (ARPES) and magnetic quantum oscil- 
lation experiments such as de Haas-van Alphen (dHvA) 
or Shubnikov-de Hass experiments identify the région 
of momentum space where zéro energy fermionic exci- 
tations exist at finite température. ARPES experiment 
directly observes the FS in the momentum space. But, 
high resolution is required to détermine the FS size. 
The quantum oscillation experiments measure the pré- 
cise value of the FS area in a spécifie plane by probing 
the oscillation frequencies of magnetization as a function 
of the applied magnetic field. The quantum oscillation 



frequency(F), the so-called dHvA frequency, is propor- 
tional to the extremal cross-sectional area Sp of the FS 
(F = hSp f'^T^e)- The quantum oscillation experiments 
also provide the information on the cyclotron effective 
électron mass m*{= {fi^ /2Tr)dSF/doj) and the geometry 
of the FS's. 

The band structure calculation is a complementary 
tool to the quantum oscillation experiment to analyze 
the complicated FS of the multiple band System. Quan- 
tum oscillation frequencies of heavy fermion materials, 
such as CeCug, UPta, and Ce(Ir,Co)In5, are explained 
well by conventional band calculations because the itin- 
érant 4/ électrons behave as conduction électrons near 
Ep- Although the geometry and volume of FS's are 
wcU explained by the DFT band calculation, the detected 
cyclotron mass m* is much larger than the correspond- 
ing DFT band mass rrn,^— because the DFT calcula- 
tion can not describe the correlated 4/ electronic states 
correctly. When the 4/ électrons are localized in the 
antiferromagnetic (AFM) compounds, such as CeRhIns, 
Celua, CeRh2Si2, the 4/- localized band model is more 
applicable to the description of the quantum oscilla- 
tion experiments j^i^ The 4/-localized band model can be 
performed by treating the 4/ électrons as core within 
the DFT (open-core DFT) band calculation,^ or by em- 
ploying the DFT-I-J7 band method [U: on-site Coulomb 
interaction).^ 

CeTIus (T — Co, Rh, and Ir) has been a prototypical 
System to study the crossover behavior between the itin- 
érant and localized 4/ électrons. CeCoIus* and Celrlng^ 
have itinérant 4/ électrons and superconducting ground 
states at low T. On the other hand, CeRhlus has lo- 
calized 4/ électrons and the AFM ground state. The 
measured dHvA frequency of each compound identifies 
the nature of Ce 4/ électrons, whether they are itinér- 
ant or localized. CeCoIng^iii^ and Celrlng^i^ have en- 
larged électron FS's due to the contribution of itinér- 
ant 4/ électrons, while CeRhlus has similar geometry of 
FS's but smaller size of FS's,^!^ For CeRhIn5, pressure- 
induced superconductivity was observed for P > 1.63 
GPa,^° and the drastic change in the FS was detected 
at a critical pressure of Pc — 2.35 GPa.^^ On the other 
hand, CeRhi_a;Coa;In5 shows the the doping-dependent 
reconstruction of FS deep inside the magnetically ordered 
state, away from the quantum critical point. The T- 
dependent évolution between itinérant and localized élec- 
trons also has been described by the phenomenological 
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two-fluid model, where the universal scaling behavior can 
be applied to various physical properties of the heavy 
fermion compounds.— 

In this Letter, we have addressed the T-dependent 
crossover from locahzed to itinérant 4/ électrons in 
Ceirins, and investigated its efFects on the FS proper- 
ties and clcctrical resistivity. The charge self-consistent 
version of DFT+DMFT approach,^^ as implemented in 
Réf. 17, is based on the fuh-potential linearized aug- 
mented plane-wave (FP-LAPW) band method.i- The 
correlated 4/ électrons are treated dynamically by the 
DMFT local self-energy, while ail other delocalized spd 
électrons are treated on the DFT level. The local 
self-energy matrix S(lli) is calculated from the corre- 
sponding impurity problem, in which fuU atomic inter- 
action matrix is taken into account,"'^^ To solve the im- 
purity problem, we use both the vertex corrected one- 
crossing approximation^^ and the continuons time quan- 
tum Monte-Carlo method i^°i^^ 

The main différence between low and high T spectral 
functions in the DFT-I-DMFT calculation is the existence 
of 4/ bands near Ep, as shown in Supplementary Fig. 1. 
Ce 4/ bands at high T are absent near Ep, and their 
spectral weights are distributed into the lower and up- 
per Hubbard bands. The spectral function near Ep can 
be well described by the quasiparticle band structures 
of other spd électrons although there is a small scatter- 
ing rate due to the hybridization between the conduction 
électrons at Ep and the incohérent Ce 4/ électrons in the 
Hubbard bands. As decreasing T, the spectral weight of 
the renormalized Ce 4/ bands is increased continuously 
(see Supplementary Movie 1.). The hybridization of the 
4/ and other spd bands produces very massive almost 
flat quasiparticle band structures near Ep. Thèse flat 
bands émerge as the narrow Kondo peak at Ep in the 
photoemission spectra.^^ 

The T-dependent FS has been extracted from the 
quasiparticle band structures. At low T, the FS's of the 
DFT-I-DMFT calculation are very similar to those of the 
DFT calculation, as shown in Fig. 1. Upon heating, the 
Ce 4/ électrons become localized and their contribution 
to Ep is suppressed. Accordingly, the areas of électron 
FS's {ai and Pi) are continuously decreased. In contrast, 
there occur rather big changes in other FS's areas upon 
heating. The areas of the g and h hole FS's on the z = 
plane grow and merge into one closed électron FS. The 
a électron FS identified at T = 10 K on the z = ir plane 
is divided at high T, and so new hole FS's appear near 
Z and R symmetry points. The continuons T-dependent 
évolution of the FS is provided in Supplementary Movie 
2. By integrating the volume of électron FS's, the oc- 
cupancy of the conduction électrons has been counted. 
It shows the continuons change from 3 to 4 électrons as 
température is decreased, which reflects the participation 
of one Ce 4/ électron in the bonding. 

Because the area of the FS is directly related to the 
size of the 4/ électron contribution to Ep, we have in- 
vestigated the T-dependent dHvA frequencies, as shown 
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FIG. 1: The T-dependent FS évolution in the 
DFT+DMFT calculation. The FS's are extracted from 
the DFT+DMFT quasiparticle band structures at 10 K(b,f), 
20 K(c,g) and 300 K(d,h). For comparison, the FS's obtained 
from the DFT band calculation are also provided (a,e). Be- 
cause the main FS's in Ceirins are nearly cylindrical due to 
the quasi-2D nature of its crystal structure, the FS's only on 
the z = and z = n planes are shown. The FS's on thèse 
planes are identified from the dHvA frequencies because the 
symmetric plane provides the extremal cross-section of the 
FS. There are two main cylindrical électron FS's represented 
by Ui and Pi branches observed in the dHvA experiment (see 
(a) and (e)). Those branches are identified at ail température 
range. On the other hand, the FS's denoted as g, h (hole 
FS's) and a, c (électron FS's) in the DFT calculation mani- 
fest topological changes with varying T in the DFT+DMFT 
calculation. Note that g, h, a, c branches were not identi- 
fied clearly in the dHvA experiments. Color represents the 
difi'erent band index. 



in Fig. 2(a). At high T, the dHvA frequencies are well 
consistent with those from the Ce 4/ open-core DFT cal- 
culation. With decreasing T, they show the continuons 
increase with the participation of 4/ électrons to Ep and 
follow the scaling behavior of ln(To/T), as shown in Fig. 
2(c). Ail the branches show the same characteristic tem- 
pérature Tq ^ 130 K. This behavior is consistent with 
the increase of the number of conduction électrons with 
decreasing T. 

The cyclotron mass corresponds to the effective mass 
of the carriers at the spécifie FS. As shown in Fig. 2(b) 
and (d), the calculated cyclotron masses also increase 
upon cooling, and follow a similar scaling behavior of ^ 
ln(To/T) with TJ" ~ 50 K. The cyclotron masses are also 
well fitted by the description of two fluid model by Yang 
et al}"^: (1 - T/To)3/2[l + ln(To/T)] with same Tg™ - 50 
K. Interestingly, TJ" is coïncident with the cohérent tem- 
pérature T* of Ce 4/ stateSf2^ but clearly différent from 
Tq . This feature reveals the 4/ électrons start to partici- 
pate in bonding through the hybridization with spd élec- 
trons at the température scale Tq , which is higher than 
the température T™ at which the cohérent heavy fermion 
electronic states are formed. Thèse results are reminis- 
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FIG. 2: The T-dependent dHvA frequencies (F) and cyclotron effective masses {m*). The dHvA frequencies (a) 
and effective masses (b) of ai and Pi branches are obtained from the DFT+DMFT method, and compared with those from 
the DFT method and experiments (exp). The corresponding FS's for each branch are provided in Fig. l(a,e), except ai that 
corresponds to the maximum frequency among Oi branches and is located between the z = and z = n planes. At high T, 
the cyclotron masses are very small, ranging from 0.4 to 0.7 mo {mo: bare électron mass) for Oi and j3i branches. Such small 
cyclotron masses are also reproduced in the 4/ open-core DFT calculation, in which only dispersive spd bands are crossing Ef- 
The low T dlïvA frequencies from the DFT+DMFT method are consistent with the results of DFT method in which the 4/ 
électrons are considered as itinérant type, (c) The renormalized AFi of each branch shows the scaling behavior of ln(ro/T) 
with the characteristic T of ~ 130 K. (d) AU the renormalized Am*'s show the similar scaling behavior, but with ~ 50 
K. 



cent of récent experiment, which shows the occurrence of 
FS reconstruction much earlier than the quantum criti- 
cal transition!^ Note that the above scaling law is con- 
sistent with the two-fiuid model^^ in which the cohérent 
4/ bands start to grow below T* . 

AU the calculated cyclotron masses at T = 10 K seem 
to be overestimated with respect to the expérimental 
values^ roughly by a factor of two. It is well known 
that the value of cyclotron mass has a substantial de- 
pendence on the applied magnctic field."* In the présence 
of the magnetic field, the effective mass can be reduced 
by the change of the hybridization, even without much 
change of the FS geometry^^ If one considers the high 
magnetic field in experiments, the calculated cyclotron 
masses would be consistent with expérimental values at 



low T. 

The continuons change of FS properties with T vari- 
ation is deeply related to the transport properties. Fig- 
ure 3 (a) provides the calculated resistivity for Ceirins 
as a function of T, which is compared to the expérimen- 
tal electrical resistivity. The electrical resistivity is cal- 
culated using the real part of the de conductivity (cr) — 
based on the DFT-fDMFT spectral function near Ep: 

^f.v ^ dw(ç^i:ï[A{\<.,uj)v^f'A{\<.,uj)v^'']. 

Here /x and v represent spatial coordinates. V , /(w), and 
V are the primitive volume, the Fermi Dirac distribution 
function, and the velocity, respectively. The calculated 
resistivities from low to high T are in good agreement 
with the expérimental resistivity. At high T, the elec- 
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FIG. 3: The magnetic part (4/ électron contribution) 
of the resistivity as a function of T. The expérimental 
electrical resistivity is obtained by subtracting tlie resistivity 
of Lalrins from tliat of Ceirins.— Inset Fig. (a), (b), and (c) 
sliow tlie broadening clianges of spectral weights at Ef at low 
(lOK), crossover (50 K), and high température (1000 K), re- 
spectively. E means the direction from M to F in momentum 
space. 

tronic carriers from dispersive spd bands become more 
and more decoupled from localized électrons in the 4/ 
shell, hence the carriers are scattered less at very high T. 
Upon cooling, the hybridization among local moments 
and spd carriers increases while the 4/ électrons remain 
very incohérent above 50 K, causing enhanced scattering 
mechanism for electric carriers. Below the scale T™, the 
électrons in the 4/ shell also gain cohérence which sub- 
stantially suppresses resistivity. Therefore, the maximum 
resistivity is observed near 50 K. Inset Fig. 3(a), (b), and 
(c) show the broadening of spectral weight at Ep, calcu- 
lated at low, crossover, and high T, respectively. The 
broadening corresponds to the scattering rate at the spé- 
cifie k-point. It is noted that the spreading of the spectral 
weight at crossover is wider than that at high or low T. 
This finding confirms that the DFT-I-DMFT calculation 
describes well the crossover behavior of Ce 4/ électrons 
with one 130 K) for the participation of 4/ élec- 

trons in the conduction and another T™ (^ 50 K) for the 



formation of cohérent heavy électron 4/ bands. 

We have examined the évolution of the heavy fermion 
state using electronic structure mcthods. As in the 
two fluid phenomenologyji^ the expérimental studies of 
other heavy fermion systems2i as well as the slave boson 
studieS j^^i^^ the crossover from the high T régime, where 
moments and quasiparticles coexist, to the low T Fermi 
liquid heavy fermion state, has a rich structure charac- 
terized by multiple energy scales. We have found that it 
is characterized by multiple scales which have a clear cor- 
respondence with physical observables. Tq is the onset 
of the sharp crossover where the small FS begins distort- 
ing towards the low T FS. At a lower T™, composite 
quasiparticles formed from /—moments and conduction 
électrons émerge, and this is signaled by a maximum of 
the resistivity. By that point, the FS has reached a shape 
which is doser to its zéro température final value, but the 
material is not yet a Fermi liquid, which is only reached at 
a much lower température Tp^. We can only put bounds 
for this quantity as being lower than 10 K for the 115 
material. 

The theory can be tested using several techniques such 
as ARPES, Compton scattering and scanning tunnelling 
microscopy, which have been developed as powerful tools 
for exploring the évolution of the electronic structure 
and are currently under wayi^ Our theory predicts that 
both T™ and increase as a function of pressure in the 
Celrlus material. More generally, it would be interesting 
to foUow thèse scales as a function of control parameters 
such as pressure and composition, to investigate the be- 
havior of T™ and in related materials which can be 
driven to a quantum critical point. 
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I. SPECTRAL FUNCTION iA(k,uj)) 

The main différence between the conventional DFT 
band structure and the Ce 4/-locahzed (open-core) DFT 
band structure is the existence of 4/ bands near Ep, as 
shown in Supplementary Fig. 1 (blue hnes). In the DFT 
band structure in Supplementary Fig. l(a), one can no- 
tice the spin-orbit (SO) split j = 5/2 and 7/2 states of Ce 
4/ flat bands, which are located around 0.3 and 0.6 eV, 
respectively (clearly distinguishable at M and A). They 
are hybridized with other conduction (spd) bands to form 
dispersive bands near Ep. Three bands containing Ce 4/ 
(j — 5/2) states are crossing Ep, and the occupancy of 
thèse bands is counted to four électrons: three spd élec- 
trons and one Ce 4/ électron. On the other hand, in 
the 4/ open-core DFT band structure in Supplementary 
Fig. l(b), there is no contribution from Ce 4/ électrons 
and only dispersive spd bands are observed. The occu- 
pancy of thèse bands is counted to three électrons due to 
the absence of one 4/ électron. 

We have also provided the DFT-I-DMFT spectral func- 
tions obtained at low (T = 10 K) and high (T = 300 K) 
température. It is seen in Supplementary Fig. l(b) that 
Ce 4/ bands at high T are completely removed from Ep, 
and their spectral weights are redistributed into lower 
and upper Hubbard bands at -2.5 eV and -1-3 eV (not 
shown here). The spectral function near Ep can be well 
described by the quasi-particle band structures of spd 
électrons only, although there is small scattering rate due 
to the hybridization between conduction électrons at Ep 



and incohérent Ce 4/ électrons at the Hubbard bands. 
Noteworthy is that the quasi-particle spd bands near Ep 
are very similar to those of the 4/ open-core DFT band 
calculation. 

At low T in Supplementary Fig. 1 (a), the 
DFT-I-DMFT spectral functions near Ep have quasi- 
particle band structures of Ce 4/ (j=5/2) states that 
are hybridized with spd states. The bandwidth of Ce 4/ 
(j=5/2) states is about 3 meV, which is ^ 1/100 smaller 
than that obtained in the DFT band calculation. That 
is why they are shown as straight lines in this narrow en- 
ergy range. The renormalized Ce 4/ bands are hybridized 
with other spd bands to produce flat quasi-particle band 
structures, as shown in Supplementary Fig. l(a). Thèse 
flat bands émerge as the narrow Kondo peak at Ep in 
the photoemission spectra. Because Ce 4/ (j=5/2) states 
are simply renormalized with respect to Ep, the quasi- 
particle bands reproduce the FS's of the DFT band cal- 
culation. 



II. THE CONTINUOUS T-DEPENDENT 
IMAGES OF THE SPECTRAL FUNCTIONS 

The images for moving pictures (Movies 1 and 2) are 
provided for the additional Supplementary information. 
Movie 1 shows the T-dependent variation of A{k,u!) on 
the k-points between F and X from T = 10 to 1000 K. 
The Movie 2 demonstrates the T-dependent variation of 
FS's on the z = plane from T = 10 to 1000 K. 
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(a) DFKbIude lines) vs. DFT+DMFT at 10 K 
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(b) open-core DFT(blude lines) vs. DFT+DMFT at 300 K 
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FIG. 1: The DFT band structures and DFT+DMFT 
spectral functions at high and low T. The 

DFT+DMFT spectral functions calculated at T = 10 K 
(a) and 300 K (b) along the symmetry lines are provided. 
The spectral fonction A{k, uj) is obtained from A{k, uj) = 
— ^lmG{k,Lo). In (a), the DFT band structures, in which 
Ce 4/ électrons are treated as valence électrons (delocalized 
limit), are drawn with blue lines. In (b), the 4/ open-core 
DFT band structures, in which Ce 4/ électrons are treated 
as core électrons (localized limit), are drawn with blue lines. 
The green lines along M and F indicate the k-point paths used 
in Inset Fig. 3 (a), (b), and (c) of the Letter. 



